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ABSTRACT 



Context. The nearby (z=0.031) TeV blazar Mrk421 was reported to be in a high state of flux activity since November, 2009. 

Aims. To investigate possible changes in the physical parameters of Mrk421 during its high state of activity using multiwavelength 

data. 

Methods. We have observed this source in bright state using High Altitude GAmma Ray (HAGAR) telescope array at energies above 

250 GeV during February 13 - 19, 2010. Optical, X-ray and y-ray archival data are also used to obtain the SEDs and light curves. 

Results. Mrk421 was found to undergo one of its brightest flaring episodes on February 17, 2010 by various observations in X-rays 

and y-rays. HAGAR observations during February 13- 19, 2010 at the energies above 250 GeV show an enhancement in the flux 

level, with a maximum flux of ~7 Crab units being detected on February 17, 2010. We present the spectral energy distributions during 

this flaring episode and investigate the correlation of the variability in X-ray and y-ray bands. 

Conclusions. Our multiwavelength study suggests that the flare detected during February 16 and 17, 2010 could arise due to a passing 

shock in the jet. 

Key words. Astroparticle physics -BL lacertae objects: individual (Mrk421)- Telescopes -Gamma rays: galaxies 



1. Introduction 

Blazars are a subclass of Active Galactic Nuclei (AGN) char- 
acterized by nonthermal emission extending from radio to high 
energies. The broadband radiation originates within a relativis- 
tic jet that is oriented very close to the line of sight. Spectral 
energy distribution (SED) of high energy peaked TeV blazars 
show two broad peaks. The first peak is located between infrared 
to X-ray energies and the second peak at y-ray energies. It is 
believed that the first peak of the SED originates due to syn- 
chrotron radiation by relativistic electrons gyrating in the mag- 
netic field of the jet. The origin of the high energy GeV/TeV peak 
is still under debate. This high energy peak might originate ei- 
ther due to interaction of electrons with photon field via Inverse 
Compton (IC) scattering as in leptonic models or due to interac- 
tions of protons with matter, magnetic field or photon fields in 
hadronic models. The seed photon field for IC scattering in lep- 
tonic models may come from synchrotron emission by the same 
population of electrons which produce the low ene rgy bump (for 
a rece nt review of observations and models, see Krawczvnski) 
( |2004|) ) in synchrotron self-Compton (SSC) model and from 
the ex ternal photons from accretion disk (Dermer & Schlickeise 
Il993h or broad line region (BLR) (GhiselUni & Madaa.igg" 
in external Compton models. The hadronic models suggest that 
high energy peak could be due to extremely energetic pro- 
tons gyrat ing in a strong m agnetic field emitti ng synchrotron 
radiation dAharon ian 2000; iMiicke et al.l l2003l) . or as inverse 
Com pton and syn chrotron emission from a proton-induced cas- 
cade (lMannheiml[l 998). 

The blazar, Mrk421 (z - 0.031), is the first extragalac- 
tic source to be detected at y-ray energy E > 500 GeV 



(iPunch et al.l 119921) . Since the first detection by the Whipple 
Observatory y-ray telescope in 1992, Mrk421 has been detected 
by various other Atmospher ic Cherenkov Te l esco pes (ACT) 
and air shower experiments ( Petrv et al. 1996t Zweerink et aj] 
1997^,'Aharonian et al.'1999';'Bose et al. '2007V Yadav et al. 200l 
Aleksic et al. 2010; Di Sciascio & ARGO-YBJ Collaboration 
120101) . and its y-ray flux has been found to be highly variable. 

Mrk42 1 was reported to be in a high state of activity during 
November 2009 to April 2010, with flaring behavior in X-ray 
(Isobe et al. 2010) and y-ray bands detected in February, 2010. 
One of the brightest flaring episodes of this source was observed 
by various experiments on February 17, 2010. Preliminary re- 
sults from different experiments show rapid flux variation in 
the very high energy (VHE) y-rays from minute to hour time 
scales. During this blight outburst, the maxi mum VHE y-ray 
flux (> 100 GeV) reached above 9 Crab units (IOngll2010l) . with 
an average flux of 4 Crab units. 

Using the newly commissioned High Altitude GAmma Ray 
(HAGAR) telescope system, we observed Mrk421 in its high 
state of activity during February to April, 2010 and also detected 
a very bright flare above 250 GeV. 

In this paper, we study the multiwavelength behavior of 
Mrk421 during its high state of activity in February, 2010 
and follow the evolution of its SEDs over a period of seven 
days, based on data from HAGAR, Fenni-LAT, RXTE-PCA and 
Swift-XRT. We also present a brief introduction to the HAGAR 
telescope and the data analysis techniques used to detect point 
like y-ray sources with HAGAR. 
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2. HAGAR 

HAGAR, an array of ACTs using wave front sampling tech- 
nique, is located at the Indian Astronomical Observatory (lAO), 
Hanle (32° 46' 46" N, 78° 58' 35" E), in the Ladakh region 
of India, at an altitude of 4270 m. The main motivation be- 
hind setting up the y-ray array at a high altitude is to exploit 
the higher Chere nkov photon density and thus achieve lower 
energy threshold (ICowsik et al.ll200lh . HAGAR consists of an 
array of seven telescopes in the form of a hexagon, with one 
telescope at the center. All the seven telescopes have seven para- 
axially mounted front coated parabolic miiTors of diameter 0.9 
m, with a UV sensitive photo-tube at the focus of individual 
mirrors. Each telescope is separated by 50 m distance from its 
neighboring telescope. The seven PMT pulses of a telescope 
are linearly added to form a telescope output called the Royal 
sum pulse. A coincidence of at least 4 Royal sum pulses out of 
seven, above a predetemiined threshold, is taken within a time 
window of 150 ns or 300 ns depending on the zenith angle of 
pointing direction to generate a trigger for initiating data record- 
ing. The HAGAR Data Acquisition (DAQ) system is CAMAC 
based. Data recorded for each event consist of 
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Fig. 1: Typical Space Angle Distribution plot obtained by 
HAGAR for a seven fold trigger. 



a. Relative arrival time of the Cherenkov shower front at each 
mirror, as measured by TDCs with a resolution of 250 ps. 

b. A Real Time Clock (RTC) module synchronized with GPS is 
used to record the absolute arrival time of these events accu- 
rate UptO jUS. 

c. The density of Cherenkov photons at each telescope is mea- 
sured by the total charge present in PMT pulses (this is 
recorded by using 12 bit QDCs). 

d. Information on Royal sum pulses are also recorded like the 
individual PMT pulses. 

e. Latch information to indicate the triggered telescopes and 
other house keeping information on various scalar readings . 

In addition, a parallel DAQ using commercial waveform dig- 
itizers with a sampling rate of 1 GS/s (ACQIRIS make model 
DC271A) is also used. The direction and energy of y-rays are 
estimated by measuring relative time delays and densities of 
Cherenkov photons at each telescope respectively. 

The performance of the HAGAR array has been studied 
by simulations, which are done in two steps: (1) Cherenkov 
emission due to y-ray and cosmic ray induced air showers in 
the atmosphere, by using the Monte Carlo simul ation pack- 
age C ORSIKA, developed by the KASCADE group JHecketalJ 
Il998i) . (2) study of the response of the array towards the 
Cherenkov radiation produced by the simulated showers. The 
performance parameters such as energy threshold, collection 
area and sensitivity of the experiment are obtained by a detec- 
tor simulation package indigenously developed by the HAGAR 
collaboration. Energy threshold of the HAGAR telescope is es- 
timated as 204 GeV in case of vertically incident y-ray showers 
for a > 4-Fold trigger condition, for which the corresponding 
collection area is 3.2 x 10^ cm^. HAGAR sensitivity is such that 
it will detect a Crab nebula lik e source at a sign ificance level of 
5cr in 15 hours of observation ("Saha e t aDl201lh . 

Cherenkov emission due to induced air showers forms a 
spherical wavefront with a large radius of curvature and thick- 
ness of ~ 1 m at the observation level. This Cherenkov emission 
mainly originates at the shower maximum region which is at a 
height of about 5 km above the ground level at Hanle. This spher- 
ical wavefront is approximated as a plane wavefront in the data 
analysis procedure, which is a good approximation at the obser- 
vation level. The arrival direction of each shower in Cherenkov 



light pool is computed by measuring the relative arrival times 
of shower front at different telescopes. The normal to this plane 
front gives the arrival direction of the incident shower. The angle 
between the direction of the shower axis and the pointing direc- 
tion of the telescope is defined as Space Angle (i/r). The Space 
Angle is estimated for every event by measuring the relative ar- 
rival time of the shower front at each telescope. 

The obser-vations were carried out by pointing all seven tele- 
scopes towards the source or background direction at a time. 
Each source run was followed (or preceded) by a background run 
with the same exposure time (typically 40 minutes) and covered 
the same zenith angle range as that of the source to ensure that 
observations were cariied out at almost the same energy thresh- 
old. Data selection was done by using parameters which charac- 
terize good quality data, in order to reduce systematic errors. 

Extraction of the y-ray signal was cariied out by compar- 
ing the ON source and OFF source space angle distributions ob- 
tained during the same night. The excess events were computed 
from 0° to Lower Limit (LL), where LL is defined as the foot 
of distant half maximum point computed by fitting a Gaussian 
function to the space angle distribution. Background space angle 
distribution was normalized with source space angle distribution 
by comparing the tails of the distributions (LL to 6.5°), since 
no y-ray events are expected in this region. This normalization 
is required to compensate for possible changes in observation 
conditions or sky conditions during ON source and OFF source 
regions. Thus the number of y-ray events are estimated as 



Number of y - rays - I.q (Si - CkBj) 



(1) 



where S ,■ and Bj are source and background events respectively 
and Ck is the normalization constant obtained using the tail of the 
space angle distribution. 

A typical space angle plot obtained for seven fold trigger 
is shown in Figure [T] The dotted (black) line histogram cor- 
responds to source observations, dashed (green) line histogram 
corresponds to background observations. The histogram shown 
as solid (blue) line corresponds to the estimation of y-ray events 
computed as the excess of ON source events over the normalized 
OFF source background events. 
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3. Multiwavelength observations and Analysis 

Data obtained on Mrk421 by HAGAR during its high state in 
2010 are combined with archival data from Ferm/-LAT, RXTE- 
ASM, Swift-BAT, Swift-XRT, RXTE-PCA, SPOL and OVRO 
for a multiwavelength study. Details of the observations and 
analysis procedure used for each data set are discussed in this 
section. 



3.1. HAGAR 

The HAGAR observations of Mrk421 were carried out for three 
months on moonless nights, during February to April, 2010 at a 
mean zenith angle of 6°. A total of 30 hours of data were col- 
lected during these obsei-vations. Suitable cuts were applied to 
the difference in the trigger rates of ON source and OFF source 
regions to ensure that observations were carried out at almost 
the same energy threshold. Some data were also rejected due to 
unavailability of all seven telescopes at the time of observations. 
After imposing such data quality cuts, a total of 18 hours of ob- 
servation is used for further analysis. Data were analyzed ac - 
cording to the procedure discussed in §2 and lBose et alJ (l2007h . 
We have analyzed only events with signals in at least 5 tele- 
scopes (> 5-Fold) in order to reduce systematic errors in our 
data. The > 5-Fold events correspond to an energy threshold of 
about 250 GeV. 



3.2. Fermi-LAT 

The Large Area Tele scope (LAT) is a pair production telescope 
(lAtwood et al.l2009h on board Fermi spacecraft. It covers the en- 
ergy range from 20 MeV to more than 300 GeV with a field of 
view > 2.5 sr. Ferm/-LAT datiQ of Mrk42 1 obtained during 2010 
February 12-22 above 100 MeV were analyzed using the stan- 
dard analysis procedure (ScienceTools) provided by the Fermi- 
LAT collaboration. A circular region of 10° radius "region of 
interest (ROI)" was chosen around Mrk42 1 for event reconstruc- 
tion from the so-called "diffuse" event class data which has the 
maximum probability of being the source photons. We only re- 
tained events having a zenith angle <105° to avoid the back- 
ground from Earth albedo. The spectral analysis on the resulting 
data set was carried out by including galactic diffuse emission 
component model (gll_iem_v02.fit) and an isotropic background 
component model (isotropicJem_v02) with post-launch instru- 
mental response function P6V3 DIFFUSE, b y using unbinne d 
maximum likelihood analysis (C ash.. 1979; .Mattox et alj|1996t) . 
A power law spectrum was used to model the source spectrum 
above 100 MeV, with integral flux and photon index as free pa- 
rameters. The flux, spectrum and source location are determined 
by using unbinned GTLIKE algorithm. 

3.3. X-ray data from RXTE and Swift 

Proportional Counters Array (PCA) (iBradtet all ll 9931) is an 
array of five identical xenon filled proportional counter units 
(PCUs). The PCUs cover energy range from 2-60 keV with a 
total collecting area of 6500 cm^. The archival X-ray data from 
PCA on board RXTE during 2010 February 13-19 was ana- 
lyzed to obtain the X-ray spectrum and light curve. We have an- 
alyzed standard 2 PCA data which has a time resolution of 16 
seconds with energy information in 128 channels. Data analy- 
sis was performed using FTOOLS (version 5.3.1) distributed as 



a part of HEASOFT (version 5.3). Data were filtered using the 
standard procedure given in the RXTE Cook Book 4 for each of 
the observations. The background models were generated with 
the tool "pcabackest", based on RXTE GOF calibration files for 
a 'bright' source (more than 40 ct/sec/PCU). The PCA spectiiim 
in the the energy range of 3 - 30 keV was fitted by using XSPEC 
with a cutoff powerlaw with line of sight absorption. The line of 
sight absorption w as fix ed to neutral hydrogen column density 
at 1.38x102" cm-2 (JDickey & Lockman 1990). 

The XRT on board Swift (Burrows et al.i i2005i) uses a graz- 
ing incidence Wolter I telescope to focus X-rays onto a CCD. 
The instrument has an effective area of 110 cm^, 23.6 arcmin 
FOV, 15 arcsec resolution (half-power diameter), and an energy 
range of 0.2- 10 keV The Windowed Timing (WT) mode data 
were used to obtain spectrum (0.3-3 keV) from Swift-XRT dur- 
ing 2010 February 13- 19. Source photons were extracted using 
a box region with the length of 40 pixels and width about 20 
pixels. Events with grades 0-2 were selected for the WT mode 
data. The spectral data were rebinned by GRPPHA 3.0.0 with 
20 photons per bin. Standard auxiliary response files and re- 
sponse matrices were used. Spectra for this source were fitted us- 
ing XSPEC version 12.3.1 with a model consisting of absorbed 
power law over the energy range of 0.3 - 10 keV 

The "Dwell" data from RXTE-ASM was obtained from 
ASM websitq3 and these data were analyzed by the method dis- 
cussed in Chitnis et al. (2009). A daily average flux between 15 - 
50 keV from Swift-BAT was obtained from BAT websitqj- 

3.4. Optical and radio data 

The optical and radio data made available in the Fermi multi- 
wavelength support program websitesQj are used in this study. 

The optical observations were made using the SPOL 
CCD Im aging/Spectropolarimeter at Steward Observatory 
(ISmith et al.ll200^ . The optical V band photometric data were 
obtained. 

The 15 Ghz radio observations were made by using a 40 

meter single dish telescope at Owens Valley Radio Observatory 

OVRO). Details of the analysis are described in .Richards et al] 



4. Results 

Mrk421 was found to be in a high state of activity during the 
entire period of HAGAR observations during February- April, 
2010, and was in its brightest state in February, 2010. The y-ray 
and X-ray fluxes decreased in the later months, but were still 
higher than that during the quiescent state. Results of observa- 
tions of Mrk42 1 using HAGAR during the three months period 
are given in Table 1 . 

Figure |2] contains the daily light curve of Mrk42 1 during 
February -April, 2010. The upper panel shows the daily aver- 
age of VHE y-ray flux obtained from HAGAR. Bottom panel 
shows daily average from ASM on board RXTE in 1 .5 - 12 keV 
It is cleariy seen in the HAGAR as wefl as RXTE-ASM fight 
curves that Mrk421 was in its brightest state in February, 2010 
in both VHE y-rays and X-rays. 

HAGAR telescope detected Mrk421 in a high state of VHE 
y-ray flux, during the period February 13-19, 2010. One of the 



http://fermi.gsfc.nasa.gov/ 



" http://xte.mit.edu/ 

^ http://heasarc.nasa.gov/docs/swift/results/transients/ 

* http://james.as.arizona.edu/ psmith/Fermi/ 

^ http://www.astro.caltech.edu/ovroblazars/index.php ?page=home 
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Fig. 2: The top panel shows the daily average light curve of 
Mrk421 during the period February - April, 2010 in VHE y-rays 
above 250 GeV from HAGAR and the bottom panel shows the 
X-ray light curve in 1.5-12 kev from RXTE-ASM, during the 
same period. 

Table 1 : HAGAR observations during high state of activity. 



Month 


Total 


Excess number 


Mean 


Signi- 


2010 


duration 


of on source 


y-ray rate 


ficance 




(min) 


events 


(/min) 


cr 


February 


479 


6418.22 


13.4±1.05 


12.7 


March 


478 


2524.23 


5.3±1.1 


4.8 


April 


116 


414 


3.5±2.1 


1.7 



brightest flaring episodes was observed on Febraary 17, 2010 
(see Figure|2]l, with the maximum flux between 6-7 Crab units. 
The source was detected with 5 sigma significance in less than 40 
minutes of observations. We investigate in the following the mul- 
tiwavelength behaviour of Mrk421 during the flaring episode. 



4.1. Flux Variability during February 10-26,2010 

The multiwavelength (radio to y-rays) quasi-simultaneous light 
curve of Mrk421 during February 10-26, 2010 based on ob- 
servations in §3, is shown in Figure |3] The top eight panels 
correspond to data from OVRO, SPOL, Swift -XRT, RXTE- 
ASM, RXTE-PCA, Swift -BAT, Fermi-LAT (0.2-2 GeV) and 
Fermi-LAil (2 - 300 GeV) respectively. The bottom panel corre- 
sponds to HAGAR data above 250 GeV. 

A clear variation of flux over a period of seven days is ob- 
served in the optical. X-rays and y-ray s during February 13-19, 
2010. The peak flux in optical. X-rays and low energy y-ray s are 
observed to be around Febraary 16, 2010 while in the VHE y-ray 
band, the peak is seen on February 17, 2010 with a possible one 
day lag compared to the situation at lower energies. The radio 
flux has not changed significantly during this time. 
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Fig. 3: Multiwavelength light cuiTe of Mrk421 during February 
10-26,2010. 



4.2. Intra-Day and spectral variability during February 
10-26,2010 

Fermi-hP^ data indicate an intra-day flux variability at ener- 
gies > 200 MeV during the TeV flare on February 17, 2010. An 
increase in the flux is seen in the first 9 hours (MJD55244.0- 
55244.4) of LAT observations, followed by a quiescence that 
lasts for a few hours (MJD55244.4- 55244.6) before the oc- 
currence of another increase in flux during MJD55244.7 to 
MJD55244.85. A similar trend was detected by the VERITAS 
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Fig. 4: Fermi-LAI and HAGAR intra-day light curve of Mrk421 
during February 17, 2010. 
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Table 2: Cutoff energy (PCA), spectral index (PCA) and photon 
index (Ferm/-LAT >100 MeV) for four different activity states 
during Febraary 13-19. 



State 



Cutoff PCA (keV) Index PCA Index LAT 



Pre Flare 

Moderate Flare 

TeV Flare 

Post Flare 



29. iH 
19.5H 
26.8H 

35.2H 



6.3 
1.5 

2.2 
2.3 



2.34±0.05 
1.91+0.03 
2.03+0.02 
2.19+0.02 



1.65+0.03 
1.67+0.15 
1.39+0.16 
2.06+0.04 



collaboration (IOn^ll2010h . HAGAR also observed a continuous 
decrease in the flux over a period of ~ 2 hrs (MJD55244.8- 
55244.86), which were simultaneous with LAT during later part 
of the night (see Figure 4). 

Spectral variation also has been detected over eleven days 
(February 12-22) of Ferm/-LAT observations and seven days 
(February 13 - 19) of RXTE-PCA observations. A photon index 
of 1.39 + 0.17, implying a flat/hard energy spectrum was ob- 
served in the Fermi-LAI data at energies above 100 MeV on 
February 17, 2010. On the other hand maximum hardening in 
RXTE-PCA data was observed on Febru ary 16 (Figure |5] and 
Table 2). S pectral hardening in th e X-rays (iTanihata et al./i2004i) 
and y-rays dZweerink et al.lll997h during strong flares has been 
reported earlier also. 



5. Discussion 

Mrk421 has shown several high states of flux activity in X-rays 
and y-ray s over the period of the past two decad es (Buckl ev et al. l 
■l996';'Gaidos et al 1996';'Krennrich et al."l999; Takahashi et al.' 



2000; Bhattetal. 2002; Blazeiowski et al. 2005; Rose et al 
2007; Giebels et al. 2007; Albert et al. 2007; Fossati et alj 
Horan et al.. ,2009; .Vemetto & Argo-Ybi Collaboration! 



2008; 



2011 



Fig. 5: RXTE-PCA light curve and spectral index are plotted in 
the first two panels, Fermi-LAT light curve in the range 0.2 - 300 
GeV and the photon index are plotted in bottom two panels. 



lBinksll201 ffl. It was observed in the brightest state of TeV y- 
ray flux during the years of 2000-2001, when i ts maximum flux 
reached as high as 14 Crab units above 1 TeV (iTluczvkont et al.l 
l201tf ). It has also shown several mild high states since 2001, but 
was not found to be in an extreme high state till 2010. Mrk421 
had brightened up again in November 2009 and was observed 
by the VERITAS collaboration in an extreme high flux state on 
February 17, 2010. The VERITAS collabora tion reporte d a max- 
imum flux of 9 Crab units above 100 GeV (OngTOlO) and the 
HESS collaboration (Tluczvkont 201 1) also found this source in 
a high state during their follow up observations. 

We have also detected this flare using the HAGAR telescope 
system, a few hours after the VERITAS observations, with a 
maximum flux of 6-7 Crab units above 250 GeV (1 Crab unit 
~ 4.2 y-rays/minute above 250 GeV). HAGAR continued the 
observations of Mrk421 over the next two months. The average 
flux during February, 2010 was found to be ~ 3 Crab units and 1 
Crab unit during March and April, 2010. 



5.1. Spectral Energy Distribution 

The Spectral Energy Distribution of TeV blazars has a typical 
two-bump shape. It is often believed that the broadband emission 
from these sources is produced by synchrotron self-Compton 
(SSC) mechanism. The lower energy bump, which peaks at in- 
frared to X-ray is interpreted as synchrotron emission emitted by 
relativistic electrons gyrating in the magnetic field of the jet, and 
the GeV/TeV bump is attributed to IC scattering of synchrotron 
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photons by the same population of electrons which produces the 
synchrotron radiation. 

A one zone homog eneous SSC model developed by 
iKrawczvnski et all (|2004|) is used to fit the multiwavelength data 
to obtain the SED. This model assumes a spherical blob of 
plasma of a comoving radius R which travels with a bulk Lorentz 
factor r towards the observer. The emission volume is filled with 
an isotropic population of electrons and a randomly oriented uni- 
form magnetic field B. Energy spectrum of the injected electrons 
in the jet frame is described by a broken power law with low- 
energy (fimin to E\,) and high-energy {E\, to fimax) components 
with indices of pi and p2. The emitted radiation is Doppler 
boosted by the Doppler factor 

5 = [r(l-/?cos(0))]-' (2) 

where /3 is the bulk velocity of the plasma in units of the 
speed of light and 6 the angle between jet axis and the line 
of sight in the observer frame. The VHE y-ray spectrum is 
corrected for absorption by the extragalactic background light 
dPranceschini et allBOOSh . 

The radius of the emission zone is constrained by the vari- 
ability time scales. Variability present in X-ray and y-ray data is 
of the order of one day, so we have chosen Ivar ~ 1 day. The 
comoving radius of the emission zone is defined as 

R-cSt,.„,/(l+z) (3) 

We have attempted to obtain SEDs for different flux states 
using multiwavelength data of the Mrk421 during HAGAR ob- 
servations of February 13 - 19, 2010 in an attempt to search for 
any changes in the physical parameters over this period. This 
multiwavelength data is divided into four states according to the 
flux state of the source as follows: State 1: Pre Flare (13-15 
February), State 2: moderate flare (16 February), State 3: TeV 
flare (17 February) and State 4: post flare (18-19 February). 
The pre flare state of the source is a moderate high state. 

The Ferm/-LAT data is divided into three bins (0.1-1 GeV, 
1-3 GeV and 3-300 GeV) to obtain spectiTim of Mrk421 for 
'State r, 'State 2' and 'State 4' by freezing the photon index 
to 1.65, 1.67 and 2.06 respectively. These photon indices were 
obtained by analyzing the 0.1- 300 GeV data from FermZ-LAT. 
The Fermi-LAT spectrum of February 17, 2010 was obtained by 
dividing Fermi-LAI data into four bins (0.1-1 GeV, 1-3 GeV, 
3-10 GeV and 10-300 GeV) by freezing the photon index to 
1.39, obtained by the analysis of the 0.1-300 GeV data. The 
best fit photon indices and cutoff energies for all four states are 
presented in Table 2. 

The magnetic field, Doppler factor (F ~ 6), electron energy 
density, break present in electron injection spectrum (E\,), pi, 
p2, minimum and maximum electron energies are used as free 
parameters while fitting the model to optical. X-ray and y-ray 
data. The angle between the jet axis and the line of sight in the 
observer frame is taken to be 2.6°. 

5.2. Evolution of the SED during fhe h/gh state 

The Fermi-LAU collaboration (lAbdo et al.l l2009l) had reported 
photon index of the energy spectrum of Mrk42 1 to be 1.78 from 
the first 5.5 months of their observations. This refers to an av- 
erage spectiTim mostly during the quiescent state. However, dur- 
ing 'State r (see Figure 6a), we find the energy spectrum to be 
flatter with a photon index of 1 .65 coiTesponding to 0.1- 300 
GeV Fermi-hP^ data. Also, we see a plateau in the light curves 
during this state in the Ferm/-LAT data (in 2 - 300 GeV energy 



bands) and 15-50 keV Swift-BAT data (see Figure 3). The 
presence of a plateau in the light curve indicates that the injec- 
tion time scales of energetic electrons are longer than the cooling 
time scales, and the entire volume of emission zone is radiating. 
Plateaux are also o bserved in other y-ray sources like 3c454.3 
(lAbdoetal.ll2011al) . 

Flares in optical. X-ray, low energy (<2 GeV) y-ray and HE 
(2-300 GeV) y-ray band were seen in 'State 2'. The coiTe- 
sponding SED is shown in Figure 6b. These low and HE y-ray 
flares could be caused by upscattering of X-ray photons by low 
energy electrons which are produced by SSC cooling in 'State 
1 ' . As an effect of low energy flare, the low energy part of SED 
in 'State 2' is much flatter than that in 'State 1'. 

The observed GeV/TeV flares above 250 GeV on 17 
February could be produced by upscattering of X-ray photons in 
Swift - BAT energy range by higher energy electrons (see Figure 
6.C for the SED). The presence of a passing shock might have ac- 
celerated the electrons to higher energies in the emission region. 

The source is brighter at lower energy y-ray in 'State 4', 
which could be due to cooling of electrons after the high energy 
flare. The corresponding SED is shown in the Figure 6d. 

The derived SED parameters for each state are listed in Table 
4, and in Figure 6e we illustrate the change that the SED under- 
goes during the four diff'erent states of the source. 

5.3. Liglit Curve 

Significant flux variation in optical. X-ray and y-ray bands is 
detected during February, 2010. The highest energy tail of the 
electron energy distribution (y > y;,^) is responsible for the pro- 
duction of the observed X-ray synchrotron continuum at > 0.5 
keV, while the TeV y-rays might be produced through upscatter- 
ing of synchrotron photons by the same population of electrons. 
The observed optical. X-ray variability during February, 2010 
may be explained by injection of fresh electron in emission zone 
and cooling of the electrons due to SSC mechanism. Swift-BAT 
light curve showed faster variability than RXTE-ASM, which 
could be due to the cooling effect of high energy electrons, which 
produce X-rays at the 15-50 keV range. 

The observed y-ray variability mainly divided into two 
bands, < 2 GeV and above 2 GeV. The 0.2-2 GeV y-rays ob- 
served by Fermi-LAI could be produced by low energy elec- 
trons through IC scattering of UV synchrotron photons. Fermi- 
LAT detected a significant variation flux in the 0.2-2 GeV band, 
observed over a period of eleven days during February 12-22, 
2010. 

The observed HE (>2 GeV) y-rays by ferm/-LAT and VHE 
y-rays by HAGAR could be produced by IC scattering of the 
electrons having a Lorentz factor in the range ~ 10'* - 10^. The 
observations from Fermi-LAI (2 - 300 GeV) show a symmetri- 
cal flare centered around February 16, 2010 but VHE y-rays flux, 
observed by HAGAR reaches peak with one day lag on February 
17, 2010. A fresh injection of high energy electrons accelerated 
by shock could be responsible for the observed flare in y-rays. 
If SSC cooling time scale is less then the light-crossing time of 
emission zone than the flare could be symmetric. 

5.4. Cross correlation study and time lag 

From the X-ray and y-ray light curves (Figures 2 and 3), we 
can see that fluxes of Mrk421 in these bands are roughly corre- 
lated. We have investigated this correlation between X-ray and 
y-ray bands using cross-correlation function (CCF) and using 
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Table 3: Correlation coefficient and time lag of Mrk421 emission during February- April, 2010 



Instruments 



Lag (Days) Correlation Algorithm Data Used 



coefficient 



2010 



RXTE- ASM vs Swift- BAT 

RXTE- ASM vs Swift- BAT 

RXTE- ASM vs HAGAR 

Swift-BATvs HAGAR 

Swift-BATvs HAGAR 

Swift-BATvs HAGAR 



0.0 


0.81 


CCF 


13-19 February 


0.0 


0.77 


ZDCF 


10-27 February 


1.3 


0.62 


ZDCF 


February to April 


1.3 


0.88 


CCF 


13-19 February 


1.3 


0.82 


ZDCF 


13-19 February 


1.3 


0.74 


ZDCF 


February to April 
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Fig. 6: SEDs of Mrk421 during Februaiy 2010. 
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Table 4: SED parameters obtained by fitting to data using t,,a,-~ 1 day 



State 


Magnetic 


Doppler 


log E„,i„" 


log £„„.,-" 


^break 


Pl 


p2 


sy,,/ 


IC„/ 


u/ 


rf 


- 


field 


factor 


[eV] 


[eV] 


[eV] 




- 


[10"] 


[1025] 


[10-'] 


[u'Jub] 


- 


(G) 


iS) 


- 


- 


- 


- 




(Hz) 


(Hz) 


(erg/cc) 


- 


State 1 


0.026 


19.5 


9.6 


12.1 


11.3 


2.4 


4.3 


1.93 


3.58 


0.9 


33.46 


State2 


0.029 


22.0 


8.0 


12.1 


11.40 


2.2 


3.9 


5.74 


7.02 


1.4 


41.83 


States 


0.029 


21.0 


9.4 


12.1 


11.45 


2.2 


4.1 


6.13 


11.53 


1.0 


29.88 


State4 


0.028 


21.0 


9.1 


12.1 


11.45 


2.3 


4.1 


5.53 


6.76 


8.5 


27.24 



*"' E,„;„: Minimum value of energy of the electrons present in the emission zone 

**' E,„ot:Maximum value of energy of the electron present in the emission zone 

*'^' Ebreak' Break in the electron injection spectrum 

*''' Sypi-iPeak value of synchrotron bump 

<'' ICpiiPeak value of IC bump 

*^' U(,:Electron energy density 

*^' 7/:Equipartition coefficient 



z-tran sformed discrete correlation function (ZCDF) (Alexander^ 
11997) a freely available FORTRAN 77 ZDCF codfl Observed 
lag and correlation coefficient between RXTE-ASM, Swift - 
BAT and HAGAR using CCF and ZDCF algorithms are listed in 
Table 3 for various cases. The time lag between X-ray and y-ray 
emissions might impose constraints on emission region and also 
it could distinguish between SSC and External Compton (EC) 
models. 



6. Conclusions 

We found satisfactory fits for all the four states with one zone 
SSC model described in §5.1. The observed pre to post flare 
evolution of SED and light curve are explained by the model. 
Changes in the physical conditions such as Doppler factor and 
magnetic field during these observations are indicated by model- 
ing of the multiwavelength data. It appears that changes in these 
parameters are related to the activity of the source. The change 
in jet flow from 'State 1' to 'State 2' may lead to a shock in 
jet. We also found small changes in particle energy density and 
magnetic field strength from 'State 1' to 'State 2'. The best fits 
obtained by fitting multiwavelength data to SSC model during 
the rising part of the flare ('State 2' and 'State 3') are found for 
power law index of the electron injection spectrum before break 
as pl = 2.2. As p 1=2. 2 is an outcome of Fermi first order mech- 
anism, it strengthens the belief that a strong shock might have ac- 
celerated the electrons in the emissi on zone and caused the flare 
by SSC mechanism. Very recently, ( Abdoet aDl2011bl) showed 
that electrons present in the emission zone of ]V[rk421 might be 
accelerated by Fermi first order mechanism at the shock front. 
The best fit parameters from the modeling of SED using mul- 
tiwavelength data for 'State 1' indicate an aged population of 
electrons with pl = 2.4 was present in the emission zone. Light 
curve during the same period suggests that the entire volume of 
emission zone was radiating. The observed change in electron 
injection spectrum pl=2.3 during the decaying part of the flare 
might be considered to be an outcome of SSC cooling mecha- 
nism. 

The observed break in electron injection spectrum could also 
be explained as a cooling break, where escape time of an elec- 
tron of Lorentz factor -/c, equals the radiative cooling time. It is 
observed from our results that the cooling break changes with 
the state of the source. We have also found that both SED peaks 
move towards higher energy as the source flux increases, and 



http://www.weizmann.ac.il/home/tal/zdcf2.html 



come back to lower energies as it decreases. Spectral harden- 
ing is also observed at X-rays and y-rays at the time of the 
flare (high state) (see Figure 5, Table 2). A departure from 
equipartition is also observed during this high state. The emis- 
sion process at the time of high activity is more complex than in 
the quiescent state. Energy dependent variability observed from 
Fermi-hP^ has shown that the source is much more complex 
than anticipated. The emission from Mrk421 also showed intra- 
day variability on 17 February 2010, in y-ray bands above 0.2 
GeV in Fermi-LAT data and >250 GeV from HAGAR data. We 
could not detect significantly this sub-hour scale variability from 
Ferm/-LAT due to the large error bars. The detection of sub- 
hour variability from Fermi-LAT and HAGAR could constrain 
the size of the emission region. The present results of corre- 
lation studies between X-ray and y-ray flux variability are not 
conclusive enough to distinguish between the one zone SSC or 
EC model. However, the observed hardening of the emission 
spectra obtained from ferm/-LAT and RXTE-PCA data with 
a one day lag (see Figure|5]and Table 2), is inconsistent with the 
SSC model and opens up the possibility of multizone SSC or EC 
models. 



7. Summary 

i. HAGAR Cherenkov telescope has detected VHE y-rays 
from the TeV blazar Mrk42 1 during 20 1 February to April, 
with the blazar being in a high state of activity above 250 
GeV. The emission reached a peak on February 17, with a 
maximum flux of ~1 Crab units, indicating a flare. 

ii. An energy dependent variation of flux in the HE y-rays 
(0. 1 - 300 GeV) has been observed by ferm/-LAT during 
the flare in February, 2010. 

iii. Multiwavelength data indicate a variation in the spectral in- 
dex during the high state. 

iv. A weak correlation between X-rays and VHE y-ray has 
been found with a time lag of ~ 1 day. 

V. The observed multiwavelength SEDs during February 13- 
19, 2010 indicate changes in the physical conditions such as 
magnetic field, Doppler factor and particle energies in the 
emission zone. The multiwavelength flare during February 
16-17, 2010 and changes in physical conditions of the 
emission zone are explained as an effect of a passing shock 
in the jet. 
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